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C.  1 


1.  Introduction 

The  solitary  wave  Is  one  of  the  interesting  two- 
dimensional  waves  that  can  occur  in  a  liquid  contained  in  a 
horizontal  channel  with  constant  rectangular  cross  section. 
It  Is  a  stable  progressing  wave  of  permanent  type  character- 
ized by  a  single  symmetrical  intumescence  whose  elevation 
decreases  to  zero  as  the  distance  from  its  crest  increases  to 
infinity.   The  amplitude  of  this  wave  is  not  necessarily  small 
compared  with  the  depth  of  the  channel.   It  was  named  by  Scott 
Russell  [1],  [2]  who  was  led  to  present  the  first  scientific 
reports  (iS^Tj  l8^4)  about  such  a  wave  after  he  observed  it 
in  the  Edinburgh  and  Glasgow  canal.   Since  then,  the  wave  has 
been  studied  in  numerous  papers.   A  brief  account  of  most  of 
the  investigations  and  a  short  history  of  the  subject  constitute 
part  of  a  paper  by  Peters  and  Stoker  [5]-   Modern  systematic 
versions  of  the  theory  based  on  the  stretching  technique  In- 
vented by  Prledrlchs  [h]    can  be  found  in  papers  by  Keller  [5]^ 
Frledrlchs  and  Hyers  [6],  Llttman  [7],    Peters  and  Stoker  [3] 
and  in  Stoker's  book  [8]. 

Heretofore  the  solitary  wave  has  always  been  regarded  as 
a  two-dimensional  phenomenon.   It  is  our  purpose  here  to  con- 
sider three-dimensional  solitary  waves  moving  in  a  horizontal 
channel  that  has  an  arbitrary  cross  section. 

In  Section  2  we  use  a  variant  of  Frledrlchs '  method  to 
secure  an  appropriate  system  of  hydrodynamlcal  equations.   This 
system  is  used  in  Section  J>   to  study  Irrotatlonal  solitary 


waves  in  a  long  channel  with  arbitrary  cross  section.   A 
potential  function  is  used  in  the  analysis  of  the  Irrotatlonal 
wave  and  consequently  the  method  of  Section  3  is  useless  for 
the  investigation  of  rotational  waves..  In  Section  4  we  return 
to  our  basic  equations  and  show  that  they  can  be  used  to  give 
an  account  of  solitary  waves  on  a  homogeneous  running  stream, 
with  arbitrary  vorticity,  in  a  horizontal  channel  with 
arbitrary  cross  section. 

In  both  the  irrotatlonal  and  rotational  cases  we  find 

that  a  first  order  approximation  to  the  wave  profile  is  given 

2  x^ 
by  the  classical  formula  q  sech  -^,    as  it  is  for  the  case  of 

a  channel  with  rectangular  section  containing  a  liquid  of 
constant  density,  without  vorticity.   In  our  cases,  however, 
as  may  be  imagined,  the  relationships  between  the  speed  of  the 
wave,  the  amplitude  q,  and  the  scale  factor  £,    are  complicated 
expressions  involving  the  solution  ^  of  a  boundary  value 
problem  that  is  of  elliptic  type.   These  relations  are  best 
left  in  the  body  of  the  text,  where  they  can  be  found  with 
other  formulas  for  the  components  of  velocity  and  pressure. 
Three-dimensional  aspects  appear  in  the  dependence  of  various 
parameters  on  the  function  f.      The  first  order  approximation 
to  the  wave  surface  exhibits  no  lateral  modification  of  ampli- 
tude; but  this  modification  does  appear  in  the  higher  order 
approximations.   The  second  order  approximation  already  pre- 
sents a  wave  surface  whose  amplitude  is  not  constant  in  a 
direction  perpendicular  to  the  axis  of  the  channel.   Thus  the 


theory  presented  here  Introduces  a  three-dimensional  solitary- 
wave  —  a  single  intumescence  with  amplitude  modulations  in 
the  direction  of  the  long  axis  of  the  channel  and  also  in  the 
horizontal  perpendicular  direction. 


2.  Formulation 

Let  a  gravitating,  incompressible,  and  inviscid  liquid, 
with  constant  density  5,  be  confined  to  a  horizontal  straight 
channel  of  infinite  length  and  constant  cross  section. 
Suppose  that  the  equilibrium  free  surface  of  the  liquid  is 
planar  and  that  it  is  parallel  i>o  the  horizontal  x,x  -plane 
of  a  cartesian  reference  frame  whose  x-,-axis  is  taken  parallel 
to  the  rigid  wall  which  forms  the  channel.   Let  the  wall  be 
defined  by  the  equation 

Xg  =  W(x  ) 

and  the  free  surface  given  by 

Xg  =  f(x^,x  ,t)  . 

The  positive  direction  of  the  Xp-axis  is  supposed  to  be  up- 
wards.  Let  g  denote  the  gravitational  acceleration,  let  p 
denote  the  pressure,  and  let  us  use  v-,  ,  Vp,  v^  to  denote  the 
velocity  components  of  a  liquid  particle.   In  terms  of  these 
quantities,  the  elementary  theory  of  hydrodynamics  predicts 
that  the  motion  of  the  liquid  is  determined  by  the  continuity 
equation 
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the  kinematic  boundary  conditions 
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plus  dynamical  conditions  at  the  free  surface;  Initial  condi- 
tions at  t  =  0;  and  conditions  at  Infinity.   The  vortlclty 
components  are  given  by 
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The  above  equations  can  be  written  In  dlmenslonless  form 
If  we  Introduce  a  typical  length  In  the  vertical  direction, 
say  h;  a  dlmenslonless  stretching  factor  /e";  and  the  dlmenslon- 
less variables 


X  =  /e  x,h~      ,      y  =  Xph"      ,      z   =  x^h"      ,      t  =  /e"  /gh  th 


U: 
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In  terms  of   these   variables   the  basic   equations  are 


(2.5) 


(2.6) 


|^=   0 
dz 

"■dx             dx        dx-" 

5u     ,     ,,    ^U          p. 

dy            dz 

=   0 


(2.7) 


V   =   w 


3^ 


(2.8) 
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The  last  equations  differ  from  those  Introduced  by  Keller 
[5]  as  an  extension  of  the  stretching  procedure  utilized  by 
Frledrlchs  [4]  to  obtain  the  two-dimensional  shallow  water 
theory  by  a  formal  perturbation  procedure.   Keller  applied  the 
stretching  process  to  the  horizontal  space  variables,  the  time, 
and  the  vertical  velocity  component.   The  equations  above  come 
from  applying  the  stretching  process  to  one  space  variable,  the 
time,  and  the  two  velocity  components  perpendicular  to  the 
direction  of  the  channel.   It  should  be  remarked  however,  that 
the  above  equations  are  In  agreement  with  Keller's  equations 
when  w  =  0  and  the  remaining  dependent  variables  are  Independent 
of  z. 

Frledrlchs  and  Hyers  [6],  Llttman  [7],  and  Lavrlentlev  [9] 
have  proved  that  the  equations  presented  above  possess  two- 
dimensional  solitary  and  cnoldal  wave  solutions  which  do  not 
depend  on  the  z -variable;  and  they  have  shown  that  It  Is 
asymptotically  valid  to  assume  for  these  wave  solutions  that 
each  dependent  quantity  can  be  developed  In  Integral  powers  of 
the  parameter  e.   For  example,  one  such  case  Is: 


h   f(x,t)=h   fl  —  , 
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00 
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In  the  sequel  we  will  be  Interested  In  the  consequences  of  the 
same  assumption  when  It  Is  applied  to  our  three-dimensional 
equations  for  motion  in  a  channel  which  is  not  required  to  have 
a  rectangular  cross  section. 


3.  Irrotational  Motion.   The  5''lltary  Wave  in  a  Channel  with 
Arbitrary  Cross  Section 

Suppose  that  the  liquid  in  the  channel  is  at  rest  until 
some  moment  when  a  disturbance  generates  a  wave  which  ulti- 
mately becomes  a  steady  wave  progressing  with  constant  velocity 
c  in  the  direction  of  the  negative  x-axis.   For  this  steady 
wave  it  is  convenient  to  use  a  coordinate  system  such  that  the 
y-axis  passes  through  a  crest  of  fixed  identity;  and  such  that 
the  mean  departure  of  the  free  surface  from  the   xz-plane  is 
zero  at  infinity.   With  respect  to  this  moving  coordinate 
system,  the  wave  form  does  not  change  with  time;  and  since  the 
liquid  possesses  zero  vorticity  when  the  motion  is  generated, 
the  motion  is  forever  irrotational.   Hence  there  exists  a 
velocity  potential  (j)(x,y,z)  such  that  for  a  liquid  particle 
with  coordinates  (x,y,z)  the  velocity  components  are 
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and  (|)  must  satisfy 
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When  u,  v,  w  in  (2.6)  are  replaced  by  (^-l)  the  momentum  equa- 
tions yield  the  energy  integral 
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where  P  is  an  arbitrary  point  in  the  liquid  and  P   is  some 
convenient  fixed  point.   Since  we  are  interested  here  in  a  wave 
of  solitary  type,  let  us  take  P  to  be  the  point  (-co, 0,0)  and 
let  us  assume  that 


(J)  (-00,0,0)  =  (j)  (-00,0,0)  =  0  ,    (t)  (-00,0,0)  = 


^ 


Tl(-00  ,Z)  =  0  , 


7r(-oo  ,0,0)  =  0 


Under  these  assumptions  [3-3)    Is  replaced  by 


(3.4) 
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an  equation  which  can  be  used  for  the  computation  of  pressure 
once  (})  is  known.   If  the  pressure  is  taken  to  be  zero  at  the 
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free  surface  y  =  t](x,z)  and  P  is  a  point  in  this  surface  then 
(3.4)  becomes 


(3.5) 
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a  condition  to  be  satisfied  at  y  =  t]  .   In  addition  to  this 
condition,  (2.8)  and  (3'1)  show  that 


(3.6)  e^^Ti^  =  (^y  -  <t)^Ti^ 


must  also  be  satisfied  at  y  =  r| .   Equation  (2.7)  becomes 


(3.7)  <l)y    -    <t>2^2    =    ° 


Which  asserts  that  the  normal  velocity  of  the  liquid  is  zero 

at  the  rigid  channel  wall  y  =  cd(z). 

Let  us  assume  next  that  (j)(x,y,z)  and  t](x,z)  can  be 
expanded  so  that 

(3.8)  ^[x,Y,z)    -  <i)Q  +  E't'-L  +  e^<}>2  +  e^'t'j  +  ••• 
and 

(3.9)  ti(x,z)  ^    er]^   +   e   t]^  +    e\^   +  ...   . 

The  idea  now  is  to  substitute  these  expansions  in  (3-2),  (3-5), 
(3.6),  (3-7);  and  then  equate  coefficients  of  like  powers  of  e. 
This  leads  to  a  sequence  of  partial  differential  systems  for 


the  determination  of  h,^^,^^,...    and  t\-.,t]^,...    .      We  need  to 
observe  here  that  the  quantity  c  /gh  in  (3-5)  must  also  be 
replaced  by  an  expansion  in  powers  of  e.   Since 
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we  have 


(3.10) 
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The  zeroth  order  system  comes  from  substituting  e  =  0  in 
(3.2),  (3.5),  (3.6)  and  (3.7).   We  find  that  ^^{x,y,z)   must 


satisfy 


i    +<t)    =  0 
oyy  ^ozz 


in  the  domain  D  (see  Pig.  1)  which  defines  the  cross  .section 
of  the  channel.   This  domain  is  bounded  by  L:   y  =  oo(z);  and 
a  segment  S  of  the  z-axis  determined  by  the  points  where  L 
intersects  y  =  0.   We  assume  that  L  is  a  piecewlse  smooth  curve 
which  crosses  y  =  0  in  just  two  points  p^  and  Pp  such  that 
t>  ^  Pp  -  P-i  ^  0.   The  equations  (3.5)  and  (3.6)  show  that 

(|)^y{x,0,z)  +(i)^^(x,0,z)  =  0 

(l)^y(x,0,z)  =  0 
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Fig.  1 
and  (5' 7)  shows  that  the  normal  derivative  i   must  be  zero 
along  L.   Since  the  above  requires  the  normal  derivative  of 
the  potential  function  ^      to  be  zero  both  along  L  and  S  we 


conclude  that 


<f)  =  F  (x) 
^o     o^  ' 


Is  the  only  solution  without  singularities  at  the  boundary. 
In  order  to  fix  Pq(x)  we  must  pass  to  the  system  of  relations 
of  next  higher  order. 

The  first  order  system  comes  from  equating  coefficients 
of  e  in  the  equations  {J>.2),    (3-5),  (3-6)  and  (3.7).   The 
result  of  doing  this  shows  that  ^-.    must  satisfy 


A    +<[,     =  -F"(x) 
^lyy   ^Izz     o   ' 


in  D.   Along  S  we  obtain  the  conditions 
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(|)-|^y(x,0,z)  =  T]-^{x,z)^^^{x,0,z) 


(5.11)     ^t^(x,0,z)    +   2 
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In  addition,  we  find  that  the  normal  derivative  i,   must  be 
'  ^In 

zero  along  L.   Ho\ 
(3-11)  shows  that 


zero  along  L.   However,  since  we  know  that  <j)  =  F  (x),  equation 


P^(x)  =  c^/ih  , 


which  means  that  we  can  take 


^o 


C  X 
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We  see,  therefore,  that  ^^    must  be  a  potential  function  in  D 
and  Its  normal  derivative  must  be  zero  on  both  L  and  S.   It 
follows  again  that  \>-.    is  a  function  of  x  alone  and,  since 
(|)-,  (-00  ,0,0)  =  c^//gh,  we  can  take 


c,x 


(j)^  =  F-^(x)  +  -i 
/gh 

where  F-]  (  -oo  )  =  0. 

The  second  order  system  shows  that  ^^   must  satisfy 

^2yy+^2zz  =  "^l^^) 
in  D  and  ^^^   must  vanish  on  L.   Along  S,  (j)  and  r]   must  satisfy 
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Thus   the   normal   derivative  of  ({>„  along   S   is  given  by 

(()„    (x,0,z)    =   -  -^  f;'(x)     . 
^2y  gh     1 


The  constant  c   can  be  fixed  by  observing  that  if  "Y   is  the 
o  '-n 

outward  normal  derivative  of  "X.  at  the  boundary  of  D  then 


(3.13  //  (X,,  +X„)dxdy=     Xn^'- 


, is 

■yy 

D  L+S 


where  s  denotes  the  arc  length  along  the  boundary  and  L+S  is 
oriented  counterclockwise  as  shoown  in  Fig.  1.   This  identity 
imposes  a  condition  for  the  existence  of  a  function  TC  which 
has  a  prescribed  normal  derivative  on  the  boundary  of  D  and 
which  satisfies 


X         +'X        =   K(z,y) 
^zz   ^^yy    ^    '-y ' 


in  D.   Such  a  function  does  not  exist  if  the  assigned  normal 
derivative  violates 


^K(z,y)dzdy  =  f     X^  ^^ 
D  L+S 


The  application  of  (3. 13)  to  (j)p  gives 
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or 
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where  Ah  Is  the  area  of  the  cross  section  of  the  channel 

occupied  by  the  liquid  at  Infinity,  and  bh  Is  the  breadth  of 
the  free  surface  at  Infinity.   Equation  (5.14)  shows  that 
either  P. (x)  =  0  or 

c^b 
(3.15)  A  -,-1^.0. 


This  Indicates  a  bifurcation  phenomenon.   The  first  alterna- 
tive. I.e.   P-,  (x)  =  0,  leads  to  a  uniform  flow  and  hence  It 
will  be  Ignored  here.   The  second  alternative  gives 

,,  -^,  2     Ah     Ah^     , 

(5.16)  c^  =  g-^=  g^^=  gd 

where  d  Is  the  mean  depth  of  the  liquid  at  Infinity.   Notice 

that  if  the  vertical  scale  factor  h  is  taken  equal  to  the  mean 

depth  then 

2  A 

c  =  gd   and   t-  =  1  . 
o  D 

This  choice  for  h  simplifies  the  equation  for  ^2'      Therefore 
we  will  hereafter  suppose  that 


h  =  d 
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In  accordance  with  the  terminology  used  for  the  case  of  a 
channel  with  rectangular  cross  section,  the  speed  c  =  /gd 
will  be  referred  to  as  the  critical  speed. 
The  function  ({)_  must  now  satisfy 

^2yy  ^2zz     1^  ' 
for  (z,,y)  in  D; 

(3.17)  (t)   (x,0,z)  =  -f1[(x) 


and 


't>o.  =  0 


2n 


along  L.   This  implies  that  we  can  write 


(3.18)  <t)p  -  -F''(x)^^(y,z)  +  F^(x)  + 
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where 


for  (z,y)  in  D; 


and 


^2zz  +  ^2yy  =  ^ 


t2y(0,z)  =  1 


^2n  =  0 


along  L.   Tliere  i?  no  less  in  generality  if  we  make  f^   unique 
cifyin.:  ■^-(0,0)  =  0.   We  assume  1 
From  eriuatjon  (3.I?)  we  see  that 


by  specifyin.:  ■^-(0,0)  =  0.   We  assume  that  ipr    is  now  known 


■  D 


(3.19) 


T1^(X,Z)  =  Tl^(x)  =  -F-J^(x) 


but  to  determine  "q,  (x)  we  must  proceed  to  the  system  of  third 
order. 

The  third  order  system  shows  that  (j)^  must  satisfy 

(3.20)      Syy^Szz  =  Pi"Mf2^y^^)-^lM 

in  D,  and  ({)-,  must  vanish  on  L.   Along  S,    the  requirements  are 
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In  terms  of  f^  ^^^   'Hn  "these  conditions  are 
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The  elimination  of  Tjp  from  (3.21)  and  (3.22)  gives 


16 


^ 
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since  the  normal  derivative  i-.  Is  zero  on  L,  the  condi- 
gn 

tlon  (3.13),  for  TC  =  ^^>    can  be  written 

P2 

D  p. 


As  we   can   see   from   (3-l8)    and    (3-23)    this   condition  requires 


D 


=   -ii'-['(x)    /      7l/^{0,z]dz  -  2hy]^ix)T]^{x) 


-Tl^(x)Tij(x)    /      ^        (0,z)dz   + 


2bc^Ti^(x) 


"gd 


-  bF2(x) 


Because  A  =  b  and  fJ(x)  =  -•n-,(x),  the  last  equation  reduces  to 


(3.26) 


r-  P' 


L  P- 


^2(0,z)dz  -  JJ   ^2{z,y)dzdy 


D 


Tin  (^) 


2bc^r||(x) 


,-     P' 


^ 


f        V'2yy(0,z)dZ+2b 


Tl^(x)ri^(x) 


L  Pi 
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which  Is  an  ordinary  non  linear  equation  for  the  determination 
of  the  first  approximation  to  the  surface  wave  t](x,z)  since  ^„ 
Is  known.   The  coefficient  of  t]'''  (x)  can  be  transformed  Into 

P-,  D 


=  /   ^2^2n^^  -o(r^2f^2yy+^2zz^^^^y 


L+S 


D 


=  If^4y^4z^^^^^ 


D 


which  shows  that  m  Is  positive.   The  coefficient  of  ii  (x)ri'(x) 


IS 


^1  -   - 


^2yy(0,Z)dZ+2b 


LPi 
p. 


^  =  - 


[1  -^2zz^°'^^^^^  +^^ 


=  -3b   1  - 


[^2^(0, p^)  -^2z^°'Pl^^ 


3b 


=  -3br 


where 


[^2z^Q-P2^  -^2z(°-Pl^^ 

-  ^  5B 
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We  remark  that  If  the  curve  L  intersects  the  z-axls  orthogonally 

then  r  =  1.   [If  ^p  (Ojp^)  -  i/^    (O^Pn  )  is  positive  one  can  make 

r  >  0  by  taking  l/b  =  d/(breadth)  sufficiently  small.]   If  in 

addition  to  the  above  definitions  of  m  and  m,  we  define  m„  as 

o      1  2 

2bc^ 

m^  =  

^    /gd 

the  equation  (5.26)  can  be  written  as 

(3.27  Vl'  ^^^  "  m^r\^{x)T]'^{x)   +m2Tij(x)  . 

This  is  the  well  known  basic  equation  for  the  theory  of  solitary 
or  cnoidal  waves  in  channels. 

For  the  solitary  wave  we  impose  the  conditions  i^,  (-od,  z)  =  0 
and  ([),  (-CD,0,z)  =  0.   The  latter  condition  guarantees  that  the 
vertical  velocity  is  zero  at  -co  ,  and  from  it,  as  we  can  see 
from  (3-17)  and  (3. 19), 

^^^{-CD  ,0,z)    =  -f!j^(-oo)  =  Tij(-oo)  =  0  . 

Also,  since  ^^    (-00, 0,0)  =  c^//gdi   it  follows  from  (3.18)  that 
FpC-oo)  =  0  which  requires  r)''(-oo),  as  (3.24)  shows.   With  these 
conditions,  namely 

Ti^(-oo)  =  Ti|(-oo)  =  r]"(-oo)  =  0  , 

and  the  symmetry  condition  t^-,  (x)  =  t],(-x)  it  is  not  difficult 
to  solve  (3.27)  and  verify  that 
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(3.28) 


Tll(x) 


Since  m  >  0,  we  have  a  solitary  wave  only  if  nip  >  0  and  this 
Is  so  only  If  c-,  >•  0. 

It  is  convenient  to  Introduce  the  quantities 


q  =  - 


5niped   mped 

m-,      "     br 


and 


£   = 


m^e 


m. 


o 


brq 

dm 
o 


where 


m 
^    A 


^(^2^4-4)dzdy 


D 


dzdy 


D 


and  we  suppose  r  >  0.   In  terms  of  these  quantities,  (3.28)  takes 
the  form 


(3.29) 


ni(x)  =-^   se 


ch 


2  xi 
2/e 


The  pressure  equation  (3-^)  can  be  written  as 


2       i 


and  by  using  the  foregoing  results  this  is  the  same  as 
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r  00 


7r{x,y,z) 


h 


k 
e   c, 


^g^ 


-  y  -  2 


OD 


kj 


i2 


e5 


r  00 


I 


k=2 


e        (f) 


ky 


+ 


oo 


£ 


k-2i 
e        (f) 


kz 


or 


Tr(x,y,z)    = 


;  EC    c^  (b 

°   +-^  -  y  -  -^  -   £(1)     (}),     +o(e) 


2gd  gd 


where  the  notation  G(e)  =  o(e)  means 


lim   V  ^  =  0  . 
e->0   ^ 


If  we   return  to   the   original   variables,    our  theory  up  to 
this   point   shows   that: 

2 

f  =  hT]   =  dT]   =  deT],    +  de  ri^  +   •  •  • 

2  -^1 
=   q   sech     -^  +   o(q) 


v^   =  u/gh  =  i^/gH  =  /^   [  <{)^^  +  e(t>^^  ]  +  o  (  £ ) 


=  /gd 


2X-1    iy 

-L   +  74:t  -  ^  seen     — ^r-r 
2d        d  2d 


+  o(q) 


v/ih        "Pj/Sd         /p-        2 
Vp  =    ^  ^     -  —^ =  y  &      r  c 

/?  /e  /£ 

3/2 


[£^<l)2y+0(£^)] 


X,  ^ 


X.  ^ 


gd    (§)  Jy  ■iJ^^{y,2.)sech     -^  tanh  -^  +  o{q^/^) 
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^3  = 


gd  (f)^^^  II  ^2z^y'^)sech2  ^  tanh  ^  +  oCq^/^) 


p  =  5gh7r  =  5gh 


2Kd 


ec  C-, 
o     o  1 


gd 


y 


ec 


•2gd 


^  [f;(x)  +-^] 
gd         /gd 


+  o(e) 


=  5g 


2  ^1-^ 
q  sech  —^  -  x^ 


+  (q) 


In  the  above  formulas:' 

q  Is  the  amplitude  of  the  first  order  wave. 

d  is  the  mean  depth  of  the  liquid  at  infinity. 


1-1 
d   d 


^  is  a  scale  factor  which  measures  the  sharpness 
of  the  crest  of  the  solitary  wave.   If  i/d  increases,  the 
curvature  at  the  crest  increases,  but  the  distance  between  the 
two  points  of  inflexion  of  •n-,{x)  decreases. 

r  and  y   are  quantities  which  depend  on  the  cross  section 
of  the  channel.   They  are  given  by 


^  =  1  -^  j   ^2zz(0'^)^^ 


ff  {i.l^^4^)aziy 


7 


D 


dzdy 


D 
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We  recall  that  V'o^y'^)  i^  defined  as  the  solution  of 
the  following  boundary  value  problem 


^2yy(y^z)  +  ^2zZ  ==  ^  '       ^^'^^      ^"^  ^ 


^2y(0,z)  =  1  ,  on  S 


^2n  "  °  '     (y.z)  on  L 


^2^0,0)  =  0 


D  is  the  domain  in  the   zy-plane  defined  by  y  =  0  and  y=(D(dz)/d 
as  shovm  in  Fig.  1.   The  rectangular  and  semicircular  cross 
sections  are  examples  of  sections  for  which  ij/r^   can  easily  be 
found.   For  the  rectangle  V'o(y>2)  =  [(y+1)  -l]/2  and  with  this, 
r  =  1  while  7  =  l/3«   For  the  semicircle  we  find 

2   2 
(3.30)  ^2^^'^^  "  ^  l^    +  y 

,2  ,.    _>2 


-  ^  Re  W^g^   in  (a+C)  -  ^^^f^   In  (a-C) 


a 

TT 


.  --2  -  4  In  a 


where  a  =  ^/ir,    C,   is  the  complex  variable  C  =  z  +  iy  and  Re  means 
the  real  part  of  the  function  in  the  brace.   From  (3-50)  we 
calculate  r  =  1  and 
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(3.31) 


7  = 


TT 


TT 


j7(^2y  +  4^^^^y 


D 
r  a 


J   V'2(0'2)dz  -  jj  ^2(y'Z)dzdy 


-a 


D 


7  = 


74  -  i^8  In  2  -  37r' 


The  velocity  of  the  solitary  wave  Is  given  by 
v-,(-aD  ,y,z)  =  c  which  Is 


qr. 


c  =  /id  [1+1^]+  o(q) 


For  the  case  of  a  channel  of  rectangular  cross  section,  this 
reduces  to  the  well  known  formula  of  Scott  Russell: 


c  =yih  (i  +  ^)  +  o(q) 


where  h  Is  the  depth. 

We  have  summarized  a  theory  which  in  some  respects 
resembles  closely  that  for  the  motion  of  a  two-dimensional 
solitary  wave  in  a  rectangular  channel.   In  fact,  if  the  domain 
D  Is  taken  to  be  a  rectangle,  all  of  the  above  formulas  agree 
with  the  known  formulas. 

Three-dimensional  aspects  are  perhaps  not  very  pronounced 
in  the  above  approximations  but  they  are  manifest  in  the 
quantities  d,  I,    Vp  and  v^  since  these  depend  on  z.   To  the 


24 


order  of  approximation  given  above,  the  amplitude  of  the  wave 
suffers  no  lateral  change.   It  is  interesting  to  note,  however, 
that  this  is  no  longer  true  in  the  higher  order  approximations. 
The  next  term  in  the  development  of  the  wave  surface  depends 
on  T]^;    and  we  can  see  from  (3.24)  that 

„         C  Tl   (x)     r]      {x) 

'  Ti  (x,z)  =  -^p(0,z)ii  (x)  +  -Li ±^ f:(x)  . 

The  term  -•)//„(  0,  z  )ti^  (x)  modifies  the  amplitude  with  respect  to 
the  lateral  z-direction.   Therefore  the  approximation 


yi^  =   f  (x^,x^)  =  deri^  \  ~d~  i  ^  ^^  ^2  \~c 


x,/e  X  \      p 


d   '  d 


furnishes  an  example  of  a  three-dimensional  solitary  wave. 

Notice  that  the  sign  of  -i/^{0,z)r]    {x)    depends  on  x.   If 
s.  and  -s.  are  the  abscissas  of  the  two  points  of  inflexion  of 

/  ^    Q     1^2  x^      q.  ^  n 

y  =  ^i(^5  =  id  ^^^^  i/i  '     ^ 

then  -q  (x)  <   0  when  |x|  <    |  s .  |  ;  and  t]'(x)  >   0  when  |x|  >    |s.|. 

For  a  semicircular  cross  section  the  lateral  change  in 
amplitude  from  the  axis  x^  =  0  to  the  side  x-,  =  k-d/ir   is 

f(x^,0)  -f(X-L  ,M)  =  d£%(0,|)Ti;[(x)+  0(£^) 


2 

¥i--^2(o4) 

^    o 


2  -^1^ 
2   sinh  -^  -  1 


4  '^1^  2 

sech  — ^  +  o(q  ) 
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where  R  =  ^d/v   is  the  radius  of  the  cross  section.   From  (5.^0) 
we  obtain 


^2(0,^) 


4(3-  4  In  2  ) 
5 

IT 


The  ratio  ■!/'o(0^:jf)/7  is 


'^2^^'   ?^  _   12(3  -4  In  2) 


7 


Therefore 


74  -  37r  -  48  In  2 


4d. 


f(x.  ,0)  -f{x   — ) 

±  -1-    TT 


24(3-4  In  2)      q^ 
7r(74-  3Tr^-  48  In  2)   ^o  >- 


2  slnh  -4--  1 
d 


4  ■^l       2 
sech  -^  +  o(q  ) 


At  x^  =0  this  is  approximately  equal  to 


f(o,o)  -f(o,  ^)  =   -(.156)  ^ 


4.  A  Rotational  Solitary  Wave  In  a  Channel  with  Arbitrary 
Cross  Section 

In  1953  Hunt  [10]  studied  the  problem  of  two-dimensional 
rotational  solitary  and  cnoidal  waves  generated  by  applying  a 
disturbance  to  an  invlscld  liquid  stream  of  constant  density 
with  vortlcity  depending  continuously  on  the  depth.   Assuming 
the  existence  of  such  waves  Hunt  showed  how  they  could  be 
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approximated.   Subsequently,  the  problem  was  Investigated  in 
more  detail  within  the  framework  of  perturbation  theory  by 
Moiseev  [11]  (I960),  Ter-Krlkorov  [12]  (1961);  and  Benjamin 
[13]  (1962).   Ter-Krikorov  proved  the  existence  of  the  two- 
dimensional  rotational  wave  by  using  the  method  which 
Friedrichs  and  Hyers  [6]  devised  for  the  irrotational  case. 
In  1959  Peters  and  Stoker  [3]  presented  an  analysis  of  a  two- 
dimensional  solitary  wave  with  a  continuous  density  stratifica- 
tion —  another  kind  of  rotational  wave.   Since  the  wave  analyzed 
in  [3]  is  assumed  to  arise  from  the  rest  position  of  the  non- 
homogeneous  liquid  its  vorticity  is  self -generated.   Shen  [l4] 
(1965)  extended  the  results  of  ['5]    and  [13]  to  two-dimensional 
rotational  waves  in  compressible  media.   None  of  the  methods  in 
these  papers  can  be  extended  to  three  dimensions  because  each 
depends  essentially  on  the  use  of  a  stream  function  peculiar  to 
two  dimensions.   In  order  to  discuss  the  possibility  of  three- 
dimensional  rotational  solitary  and  cnoidal  waves  we  must  return 
to  the  basic  equations  (2.1) -(2. 4). 

The  purpose  of  this  section  is  to  show  that  the  variant 
of  Friedrichs'  technique  which  leads  to  equations  (2.5)-(2.8), 
plus  some  auxiliary  devices,  can  be  used  to  investigate  rota- 
tional solitary  waves  in  a  stream  of  constant  density  with  an 
arbitrary  lateral  distribution  of  velocity,  v  (y,z),  (inde- 
pendent of  x)  when  the  stream  is  confined  to  a  channel  of  arbit- 
rary cross  section.   We  assume  that  a  coordinate  system  is  chosen 
as  it  was  in  Section  3;  and  we  suppose  that  the  vertical  scale 


27 


factor  Is  taken  equal  to  the  mean  depth  of  the  liquid  at 
infinity.   That  Is,  as  before, 

h  =  d  . 

With  respect  to  the  above  coordinate  system,  the  basic  equa- 
tions are: 


(^.1) 


Siu    Sv    ^w  ^  „ 
Bx    ^    ^  " 


(4.2) 


/■   r   Su  ,  StTh  ,    Su  ,    5u    ^ 


{ 


r    Sv   ,   T   ,   StTt  ,     Sv  ,     Bv 


=  0 


5w  ,  ^ 


Bw  ,    S 


^[-S-S]--^--§f  =0 


(4.3) 


V  =  w 


for 


y  =  co(z) 


(4.4) 


eu 


^=v-w^   for   y  =  Ti(x,z) 


(4.5) 


e7r(x,T],z)  =  0 


The  Idea  once  again  Is  to  develop  the  dependent  quantities  in 
integral  powers  of  e  and  then  solve  the  partial  differential 
systems  which  result  from  equating  coefficients  of  like  powers 
of  e. 

The  horizontal  velocity  component  u  —  the  axial  velocity 
is  to  be  regarded  as  a  sum  of  two  parts 
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Vq ( y , z ) 

u  =  +  U(x,y,z)  . 

/gd 

The  velocity 

Vq(x2,x^)  =  V^(dy,dz)  =  VQ(y,z)  , 

which  we  suppose  is  continuous  and  non-negative,  is  the  velocity 
of  the  undisturbed  liquid  stream  which  has  no  lateral  velocity 
components,  and  which  is  confined  to  a  horizontal  channel  of 
infinite  length  with  arbitrary  cross  section.   We  assume  that 
the  velocity  U  due  to  a  steady  state  solitary  wave  can  be 
expanded  in  the  form 

c  p 

U  =  +  eu-  +  e  Up  +  .  .  . 

/gd 

We  can  thus  suppose  that 

2 
u  =  u  +  eu^  +  e  Uo  +  •  .  • 
o    1     2 

where 

c^+  v^(y,z) 


(4.6)  u^  = 


o   o 


^ 


Guided  by  what  we  have  found  for  the  irrotational  case  we  also 
assume  that 

2     3 

V  =  e  Vo+  £  V-,  +  .  . . 

2,3^ 
w  =  e  iXp  +  e  i-L^  +  . . . 

2 

TT  =  -y  +  ETT-,  +  £  TTp  +  .  .  . 


2 


29 


If  we  substitute  the  above  expansions  In  (^.1),  (4.2), 

2 
(^•3)>  {^-^),    (^'5)  and  equate  coefficients  of  e   in  each 

equation  we  find  the  system: 


(4.7)  ^lx+^2y  +  ^2z  =  ° 


(4.8)  u  U-,  +17-,  +  v„u^  +  )j,„u^„  =  0 

^    '  o  Ix   Ix   2  oy  ^2  oz 


(4.9)  TT^y   =   0 


(4.10)  TT^^  =  0 

(4.11)  ^2  =  ^^2  If  ^   y  =  ^^^^ 

(4.12)  Uq(0,z)ti^^(x,z)  =  V2(x,0,z) 

(4.13)  7r^(x,0,z)  =  T)^(x,z)  . 

Equations  (4.9),  (4.10)  and  (4.13)  show  that 

7r^(x,y,z)  =  r]^{x) 

where  t^,  (x)  is  a  function  of  x  only.   The  elimination  of  u^ 
from  (4.7)  and  (4.8)  gives 

Let  i_  and  j_  be  unit  vectors  along  the  z-axis  and  y-axis 
respectively.   Let  V  be  the  vector 
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and  let 


V  =  111  +  vj 


'-i^-ii 


If  n  is  the  unit  outward  normal  of  the  domain  D  defined  in 
Section  J>,    the  divergence  theorem  states 


(^.15) 


11    ^'J.   dzdy  =  /    V  .  n  ds 
D  L+S 


Along  L  we  have 


{\\.^1_   +  V2J_).n  =  0 


as  we  can  see  from  (^.11) 
vector 


Hence  if  we  apply  (4,15)  to  the 


the  equation  (4.14)  implies 


P, 


'Ix 


dzdy 

w.  — 2~ 
■^   u 
D    o 


V2(x,0,z)dz 


u 


P- 


and  then  from  (^.12)  we  have  the  condition 


^ix^^^ 


dzdy 


D   ^o 


-  b 


0 


Here  again  we  have  a  bifurcation  phenomenon.   Either  t)   =  0  or 

_L-A. 

else 


D   ^o 
31 


=  b 


The  first  alternative  Is  uninteresting  since  It  allows  only 
the  equilibrium  flow.   The  second  alternative  Is  an  equation 
for  the  determination  of  c  ,  namely, 

(4.17)  ff ^^ ^  =  ^. 

D   t^o+^o^^^y^j 


Equation  (4.17)  specifies  two  values  of  c  .   This  can  be 
seen  by  tracing  the  change  In  the  Integral 

as  c   varies  from  -oo  to  oo  .   We  suppose  that  v  (z,y)  Is  a  non- 
negative  continuous  function  on  the  closure  of  D,  and  that  m 
and  M  are  respectively  Its  minimum  and  maximum  values  there. 
We  cannot  admit  values  of  c   between  -M  and  -m  because  for  such 
values  l(c  )  falls  to  exist.   If  c   varies  from  -m  to  +00,  I(c_) 
changes  from  +00  to  zero;  furthermore  as  we  can  see  from 

i'(o„)  =  -2  ff S5ay__^ 

it  never  increases  In  the  range  -m  ^  c  .   Hence  there  is  one  value 

^        o 

c  >   -m  which  satisfies  (4.17);  and  this  value  corresponds  to  a 
wave  advancing  upstream.   If  c^  varies  from  -oo  to  -M,  1(0^) 
changes  from  zero  to  +co  ;  and  in  the  range  -00  *^  c   -^  -M, 
I  (Cq)  ^  0.   Thus  there  is  one  value  of  c   <  -M  which  satisfies 
(4.17)  and  it  corresponds  to  a  wave  progressing  downstream. 
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These  values  of  c   will  be  referred  to  as  critical  values. 
More  can  be  said  about  equation  (^.17)  and  we  will  continue  the 
discussion  in  the  next  section. 

In  order  to  find  the  dependent  quantities  u,  ,  \x^,    v„  and 


T]     we  turn  to  the  next  partial  differential  system  which  comes 
from  equating  coefficient; 
tlons  (4.1)-(4.5).   It  is; 


2 

from  equating  coefficients  of  e   in  each  of  the  expanded  equa 


(^.18)  Ug^  +  v^y+H^^  =  0 


(4.19)    ^o^2x  +  ^l^lx+^2x  +  ^2^1y+^3V^^2^1z  +  ^3''oz  =  ° 


(4.20)  u^v„  +7r„   =  0 

*     '  o  2x   2y 


(4.21)  u  p.„  +  7r„  =  0 

*     '  o^2x    2z 


(4.22)  v^  =  M.^a)'(z)  ,    y  =   go(z) 


(4.23)    u^^{0,z)r]^T]-^,^  +  u^{x,0,z)T]^^+u^{0,z)T]^^ 

=  v^^{x,0,z)t]^   +v^(x,0,z) 

(^.24)  7r2(x,0,z)  -  Ti2(x,z)  . 

The  differentiation  of  (4.14)  with  respect  to  x  gives 
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^  o 

and  if  we  use  (4.20)  and  (4.21)  the  last  equation  becomes 


'^  \  U  /        \  u  /       u 

o       ^  o        o 


We  also  have  from  (4.20)  and  (4.12)  that 


(4.27)      Tr2y(x,0,z)  =  -u^V2^(x,0,  z)  =  -%(0^z)ti^^(x) 


In  addition,  (4.11)  shows  that 


or 


^2x  =  ^^2x^'(^) 


u  v^   =  u  M,„  cd'  (z) 
o  2x    o^2x  ^     ' 


^2y  "  ■^2z^'(^) 


which  means  that  the  normal  derivative  tt^  must  satisfy 


(4.28)  TT^^  =  0 


along  L. 

The  solution  of  the  boundary  value  problem  posed  by  (4.26), 
(4.27)  and  (4.28)  can  be  expressed  In  the  form 

(4-29)  7r2(x,y,z)  =  -^(  z,y  )ti^^^  +  F(x)  . 

Here  ^(z,y)  must  satisfy 
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O  0       0 


for  (y,z)  in  D.   This  is  an  elliptic  partial  differential  equa- 
tion, but  not  Laplace's  equation,  since  u  ,  as  given  by  (4.6), 
is  a  known  function  of  y  and  z.   The  boundary  conditions  for  f 
are 

(4.31)  ^y(0,z)  =  u^(0,z) 


for  the  segment  S,  while  along  L 
(4.32)  f^=0 


n 


We  will  make  ip   unique  by  imposing  (without  loss  of  generality) 
some  convenient  additional  condition  at  a  point,  say 

(4.33)  ^(0,0)  =  0  . 


The  functions  Vp,  |Xp  and  u.,  can  now  be  expressed  in  terms 
of  f   and  T)  .   From  (4.20)  and  (4.21)  we  have 


^o 


^2x  =   -^2y  =   Vlxx 


%^^2x  =  -^2z  ^  ^z^lxx 


and  integration  with  respect  to  x  gives 

%^2  =  Vlx  +  ^l^^'^^ 
^o^^2  =  ^z^lx+^2(y'^^ 
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since  we  are  Interested  here  primarily  In  the  solitary  wave  we 
impose  the  conditions 


(4.3^) 


Ti^(oo  )  =  0  ,    u^(oo,y,z)  =  c^ 
VgCoo^y^z)  =  (^^(oo^y.z)  =  0  . 


If  V2(oo,0,z)  =  0  equation  (4.12)  requires  r\      (oo)  =0  and 
consequently 


(4.35) 


^2^^'^'^^  -   ujy,z)  ^ix^^) 


(4.36) 


'>PJy,z) 

^2^^'^'^^    =   u  (y,z)  ^lx<^) 


For  the  function  u^,  equation  (4.7)  requires 


u 


Ix 


-^2y-^^2z 


[57  \u^  }      ^\u^ 


'Ix 


and  after  integration  with  respect  to  x  and  use  of  (4.^4) 


(4.37) 


u-L  -   c^ 


t 


y 


"Sy  \  u  /  'Sz 


+ 


f. 


u 


If  we  define  functions  P  and  Q  by 


(4.38) 


^^(y^z)  ^„(y,z) 


we  have 
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(^.39) 


VgCx^y^z)  =  P(y,z)Ti^^(x) 


(4.40) 


|j-2(x,y,z)  =  Q(y,z)Ti^^(x) 


(4.41)      u^(x,y,z)  =  c^  -  [Py(y,z)  +Q^(y,z)]Ti^(x)  . 

To  this  point  we  have  not  used  the  equations  (4.l8), 
(4.19),  (4.23)  and  (4.24).   We  proceed  to  show  that  these  equa- 
tions serve  for  the  determination  of  t),(x).   The  elimination 
of  Ug  from  (4.l8)  and  (4.19)  leads  to  an  equation  which,  after 
some  manipulation,  can  be  written  in  the  convenient  form 


(4.42) 


S 
^ 


3 


^2^1 


u 


o    u 


+ 


^ 


M.. 


u. 


M-g^l 


o    u 


o  -' 


v 


2x 

■~? 

u 
o 


2u,  tTt 
1  Ix 


U' 


o 


If  we  set 


V 


[u^\x^-\i^-a-^]    _   [u^v^-VgU^] 


1  +■ 


u. 


o 


u 


and  apply  the  divergence  theorem,  (4.15),  to  (4.42)  we  find 


(4.43) 


D 


IT 


2x 


2uTTr^ 
1  Ix 


u 


o 


U' 


dzdy  - 


"2 


dz 


u 


The  boundary  conditions  (4.23),  (4.24)  and  (4.12)  show  that 
(4.43)  is  the  same  as 
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(4.W)  ff 


D 


2x     1  Ix 


—2 J 

o      o 


dzdy 


[u^y(0,z)Ti-^Ti^^+  u^(0,z)7r2^(x,0,z)  -  v^y  (x,  0,  z  )ti^  ] 


Uq(0,zJ 


dz 


Then,  after  using  the  results  (4.29),  (^-39),  (^-^1)  and 
iT^    =  'n-(x),  we  have 


(4.^5) 


-^^1XXX+^'     g[C-^-(Py+Qz)^l^^lz,  .  , 
2 ^ ^  ^^^y 

D   I      o  o 


P. 


u  {0,z)ri^r]^  P  (0,z)-n-,TiT  ") 


u  (0,zj 
o 


'Ixxx 


u^(0,zj 


The  coefficient  of  f'(x)  vanishes  by  virtue  of  (4.l6)  and  there- 
fore (4.45)  reduces  to 


P. 


(4.46)   j  /   ^(0,z)dz  - 


%^  ^^^n  -iixxx 


D 


u. 


[u   (0,z)-P  (0,z)] 
u  (0,z) '^^~ 


2(P  +QJ     ) 

^^3 dzdy  n^T)^^ 


D      o 


+ 


f-.  If 

^  D 


T-f^ix 

O   y* 


dzdy 
u 


which  is  again  the  familiar  ordinary  non-linear  differential 
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equation  for  the  determination  of  the  first  order  approximation 
to  the  surface  wave  r)(x,z)  of  solitary  type. 

The  coefficient  of  r]  can  be  transformed  as  follows. 

Po 


^(0,z)dz  -  IJlL^^   dzdy 


o   J  _ 

P^  D    % 


P2 


,(0,z)dz-//,f^(-|).^(-|)|dzdy 


P^  D    ^     -o        "o 

=  J        ^(0,z)dz  -  j   ^(o,z)  ^ dz  +  yj  -^L.^-^-  dzdy 

P-L  P-,^         %(0.z)      '.D     Uq 

2 
From  (4.51),  namely  ^^y^"^'^)  ""  Uq{0,z)  we  see  that 

m^  =  //  — ^^—5 dzdy 

D      o 


which  shows  that  m  is  always  positive.   With  respect  to  the 
coefficient  of  'Hn'n-,   notice  first  that 


*>•      o         o  ^ 


+ 


oy^y   oz  z 

2 

u 
o 


1  ,   oy'^y   oz^z 
^   2 


Then 
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m,  = 


P- 


u   (0,z) 
oy  '  ' 

Uq(0,z) 


^_(0,z)    u  ^  fO,z)^,  fO,z) 


yy 

u^(0,z) 


+ 


oy 


y 


u^(0,z) 


dz 


-  2 


dzdy 


u 


o 


+ 


^,,(0,z) 

^ dz  - 


Pi   o 


u^(0,z) 


D   % 


^. 


y 


+ 


a  /*. 


o         o 


dzdy 


^2  n 


'1   L. 


2u^y(0,z)^y(0,z)    ^yy(O.z) 

^ + 


uj(0,z) 


u^(0,z)    uj(0,z) 


dz  -  3 


dzdy 
D    o 


m,  = 


a 

^ 


"^2(0, z)- 


u^(0,z) 


dz  - 


D    o 


Note  that  if  the  free  surface  of  the  liquid  at  Infinity  inter- 
sects the  channel  wall  orthogonally  then 


a 
^ 


>2(0,z)' 


u^(0,z) 


dz  =  0  . 


The  coefficient  of  t]    will  be  denoted  by  m^,  that  is. 


ITlg  =  2c^ 


dzdy 


D    o 
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In  terms  of  the  coefficients  Just  Introduced  the  equation 

(4.^6)  Is 

For  the  solitary  wave  we  Impose  the  conditions 

Tlj^(OD,z)  =  0    ,    Uj^(oo,y,z)  =  c^ 

{kA8) 

v^(oo,y,z)  =  0  ,   p,^(cD  ,y,z)  =  0   . 

As  we  have  already  seen  these  conditions  Imply  t]  (oo  )  =  t]      (oo)=0 
From  (4.24)  and  (4.29)  the  second  order  profile  r|p(x,z)  Is 

Ti2(x,z)  =  iT^{x,0,z)    =  -^(0,z)t]^^^(x)  +F(x)  . 

In  order  to   satisfy 

0  =  ri2(oo,z)    =   -tlj{0,z)r\^^co  )  +F(od  ) 

we   take  F{oo)    =   0  and  r]^      (oo  )    =   0.      Subject   to   the   conditions 

rij(oo)    =  ri|(oo  )    =  Ti^(oo  )    =    0 

and  the  symmetry  condition 

Tl3^(-X)  =  Tl^(x)   , 

the  solution  of  (4.47)  Is 

(^.^9)  Ti,(x)  =  -  -^  sech^  J  P^ 

'1   '      m,       2  Jm 
1         ''  o 


.41 


provided   c.    Is   chosen   so   that   nip  >•   0. 

For  later  use  It  Is  convenient  to  define 


and 


where 


i  = 


fm^e 


m. 


o 


K   = 


dzdy 


D 


u. 


1 
3 


o 


^m^ed 
inZ 


q 

d 


-m-, 


5% 


S 
^ 


>2(0,z) 


u  (0,z) 
o 


dz 


In  terms  of  these  quantities  (4.49)  Is 


(4.50) 


•n^  =  -i-  sech' 
'1    ed 


2/e 


We  could  go  on  to  higher  order  approximations  but  we  do 
not  think  that  this  Is  the  place  to  show  how  the  ensuing  com- 
plications can  be  dealt  with.  Instead,  It  Is  more  to  the  point 
to  collect  the  results  we  have  found  and  express  them  In  terms 
of  the  original  variables  with  respect  to  a  coordinate  system 
moving  with  the  wave.  We  have  shown  that  the  free  surface  Is 
given  by 

f  =  dT^  =  deT^,  4o(e) 

X  £, 
=   q  sech^  -^+  o(q)  . 

The  long  Itudlnal  velocity  has  the  form 
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V,  -  /p" 


v^(y,z)+c^ 


/gd 


+  eu^ 


+  o(e) 


V^{x2,x^)+c^  +  If 


Kq 


dzdy 


D   % 


-  q 


||A.^°y-yy°^-^(3eoh^^.o(q) 


i 


The  lateral  velocity  components  are  given  by 


^2  =  J 


Mv  =  fM  [e2v,  +  o(B2)] 


3/2  ,  ,^  Nl/2  ^  (y,z) 


^(  ^r"(J^r"  ^^^;;r^  sech^  ^.  tanh  ^.o(q^/2) 


m. 


u^TyTzJ 


2d 


2d 


and 


^3^ 


a\^/^  /  K  \^/2  ^  (y,z)     pxi       xi      / 
"gdl^l     —I    -^r r  sech"^  -i^  •  tanh  ^+o(q^/'^) 


2d 


2d 


The  pressure  p  Is  given  by 


p  =  6gd7r  =  6gd  [ -y  +  STr-,  +  o  (  e)  ] 


2  ^2. 
=   5g  [q  sech  -^  -  Xg]  +  o(q) 


and  hence  to  this  order  of  approximation  the  pressure  obeys  the 
hydrostatic  law. 

These  formulas  summarize  the  results  of  this  section  and 
in  fact  subsume  the  results  of  the  last  section. 
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For  ready  reference  we  append  a  table  of  definitions  of 
the  terms,  other  than  those  with  obvious  meanings,  in  the  above 
formulas: 

q  is  the  amplitude  of  the  first  order  wave. 

d  is  the  mean  depth  of  the  liquid  at  infinity. 


n  _  1  rw  . 


d   d  Jdm 


is  a  scale  factor  which  is  a  measure  of  the 
sharpness  of  the  crest  of  the  solitary  wave.   If  -r-  increases 
the  curvature  at  the  crest  increases,  but  the  distance  between 
the  two  points  of  inflexion  of  the  solitary  wave  decreases. 


dzdy   1 


D   ^o 


5 
^ 


■>P^{o,^) 


u^(0,z) 
o 


dz 


"^o  = 


D 


— ^  g  ^  dzdy 


u 


o 


TpiYfZ)    is   defined  by 


3    (%\^    SC*= 


u. 


u 


V'y(0,Z) 


1 

u 
o 


u^(0,z) 
o 


(y,z)    in  D 


^     -   0 
^n 


(y,z)    on  L 


V/(0,0)    =    0      . 

D  Is  the  domain  in  the  zy-plane  defined  by  y  =  0  and 
y  =  cd{z)  =  W(dz)/d  as  shown  in  Pig.  1. 
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/gd  /gd  /gd 

where  V^[y.,^,xJ)    is  the  longitudinal  velocity  of  the  undisturbed 
liquid  stream. 

c   Is  a  critical  speed  specified  by 

/T      dzdy      ^  Jo_ 
^^    [c^+v^(z,y)]^  ~  S^ 

where  bd  is  the  breadth  of  the  free  surface  at  infinity. 
The  velocity  c  of  the  solitary  wave  is  given  by 


-c  =  [v^(-oo ,y,z)  -  v^(y,z)]  , 


that  is. 


(4.51)  -c  =  c^  +  -r^-^—  +  o(q) 

D   ^o 


5.  Discussion 

Three  dimensional  aspects  in  the  approximations  given  in 
Section  4  are  manifest  in  the  quantities  d,  H ,    v, ,  v  and  v  ; 
but  to  the  order  of  approximation  given  above,  the  amplitude  of 
the  wave  suffers  no  lateral  change,  i.e.  no  change  in  the  z- 
dlrection.   It  is  interesting  to  observe  however,  that  the 
higher  order  approximations  do  show  changes  in  amplitude  in  the 
z-directlon.   The  next  term  in  the  development  of  the  wave  sur- 
face depends  on  t]^;    and  we  can  see  from  (4.24)  and  (4.29)  that 
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Ti2(x,z)  =  ir^{x,0,z)    =    -Tp{0,z)r]^^^{x)   +F(x)  . 

The  term  -tpiO ,z)-r]-.       (x),  whose  sign  depends  on  x,  modifies  the 
amplitude  with  respect  to  the  lateral  z-direction.   Hence  the 
approximation 

f  =  d£T]^  +de^T)2+  o(e^) 

evidently  represents  a  three-dimensional  solitary  wave. 

The  linear  theory  of  surface  waves  on  a  liquid  stream  in 
a  channel  is  based  on  the  equations  which  can  be  obtained  from 
(4.1)-(4.5)  by  substituting  e  =  1,  u  =  u^  +  u*,  v  =  v*,  w  =  w*, 
TT  =  -y+TT*;    1-1  =  11*  s.n6.   then  neglecting  products  of  starred 
quantities.   The  equations  are: 

(5.1)  u*  +  v""  +  w""  =  0 

^  -^   •'  X   y   z 


U  U*  +V*U    +  W*U  =  -IT* 

ox     oy     oz  X 

(5.2)                <   U  V*  =  -TT* 

\-j      J                          ^   o  X  y 


U  W*  =  -IV* 

ox  z 


(5.3)  V*  =  w'^a)'(z)    for   y  =  co(z) 

(5.4)  u^Ti*  =  v*      for   y  =  0 

(5.5)  -T^*  +  TT*  =  0   for   y  =  0  . 

If  we  eliminate  u*,  v*,  w*,  and  r|*  we  find  that  tt*  must 
satisfy 
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a 
^ 


^  U_  ''     "^  ^  u 


+ 


a         z 


=  0 


for  {y,z)  in  D  and  the  boundary  conditions 


UQ(0,z)7r^(x,0,z)  +  7r*(x,0,z) 


=  0 


for  S  and 


IT      -    0 
n 


for  L.   For  solutions  of  the  form 


(5.6) 


TT*  =  ■>p''{Y,z)e 


IVX 


we  must  have 


(5.7) 


o         o     o 


^*(0,z)  =  vV(0,z)^*(0,z) 


(y,z)  in  D 


^n 


{y,z)  on  L 


The  divergence  theorem  applied  to  (5-7)  yields 


(5.8) 


D 


^  dzdy 

u 
o 


T//  (0,z)dz 


As  V  — >■  0  the  wave  length  of  (5-6)  — >  oo  ;  but  f     — >■  constant, 
The  condition  (5' 8)  then  becomes 
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(5.9)  ff^^ff (gd)dzdy    ^  ^ 


which  we  can  expect  to  yield  the  approximate  speeds  of  long 
waves  progressing  against  and  with  the  stream.   Equation  (5.9) 
is  of  course  precisely  (4.17)  which  gives  the  approximate  speeds 
of  the  solitary  wave.   For  a  rectangular  channel  (5-9)  becomes 

(5.10)  /%=! 

-1  ^O 

or 

(5.11)  f   3Z ^   1. 

The  last  equation  was  derived  by  Burns  [15]  from  a  two- 
dimensional  analysis  based  on  ideas  different  from  those  above. 
Burns  gave  a  method  for  estimating  the  speed  of  long  two- 
dimensional  waves  relative  to  the  mean  of  the  velocity  v  (y) . 
His  method  is  also  applicable  to  the  equation  (5' 9)  as  we  pro- 
ceed to  show. 

If  the  mean  flow  is 


^o  =  i  ff^o^^'^^'^^'^^ 


D 
then 


v+c  -.rp(v+c)  ., 

Also   since  b/A  =   1,    (5-9)    Is   the   same  as 
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D   ^o 


=  A 


and  therefore 


(v  +  c  ) 
^  o   o^ 

gd 


1 


u^dzdy 


D 


A- 


n2 


u  dzdy 


L  D 


dzdy 
— 2~ 

D    ^o 


A- 


I  %  I dzdy 


D 


dzdy 


D   l^ol 


Successive  applications  of  Schwartz '  inequality  produces 


(v  +c  )^ 

^  o    o'   J,  1 


gd    -^3 


I Uq I dzdy 


D 


^  D  '  o 


dzdy 


-A^ 


rr     dzdy  . 

J  J     lu  |V2 

D   '  o '     1-  D 


u^l  /   dzdy 


^2 


^^ 


u^l^/^  dzdy 


D 


-D   1% 


dzdy 


D   1% 


dzdy 


UqI  ^   dzdy 


^  D 


^^ 


u^l  /   dzdy 


D 


•-  D 


u 


dzdy 
1/32 


A- 


u^l  /  P  dzdy 


D 


dzdy 


•-  D   u. 


1/2 


2p+l 
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Now 

lim   |x|^/'^  =  1  . 
n — >  00 
Hence 


(5.12)  l^o+'^ol  -^^   • 


That  Is,  the  magnitude  of  the  critical  speed  with  respect  to 
the  mean  velocity  of  the  flow  Is  never  less  than  the  critical 
velocity  of  the  flow  without  vortlclty.   The  same  thing  holds 
for  the  first  order  wave  velocity  as  we  can  see  from 


'<q  J  f 

V    -c  =   v+c^+  ^s +  o(q)    . 

°        °  ^^dzdy 


If   c      >   -m,    c    +  m  >   0  then 
o  o 


v+c      >m+c      >   Q) 
o        o  o 


c     +v     >c     +m>0 
o        o  o 


K    >   0 


.   "T 


^//^^° 


and  therefore 


If   c     <    -M: 
o 


D        ^o 


v    -  c   >  Vgd 
o  ^  ^ 


V    -  c      >   0 
o        o 


V    -  c^   >  /gd 
o        o  —  ^  ^ 
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—  —  '*^Q  JTT 

v^+  c  =  v^  -  c^ —  +  o{q) 

o      0   0     rr  ^   ^ 

o 


c^+v^  <  .M  +  v^  <  0 


.ff-^<0 


D   ^o 


K  ^  0 


and  therefore 


v^  +  c  >  ygd   . 
o      ^  ° 


If  the  amplitude  q  and  the  channel  section  D  are  given, 
the  actual  computations  of  the  wave  profile,  particle  velocity 
and  pressure  depend  on  the  evaluation  of  quantities  which 
involve  the  function  f.      This  is  a  considerable  task  and  it  is 
even  more  difficult  to  estimate  the  effects  of  vorticity  by 
comparing  parameters  in  the  rotational  wave  with  the  corre- 
sponding parameters  in  the  irrotational  wave  which  depend  on  ^p. 
We  do  no  three-dimensional  computation  here.   Instead,  we  will 
devote  a  little  attention  to  the  more  tractable  case  of  the  two- 
dimensional  rotational  wave. 

For  a  channel  with  a  rectangular  cross  section  of  depth  h 
containing  a  stream  whose  velocity  V  (xp,x-,)  =  V  (hy,hz)  =  v  (y,z) 
is  a  function  of  y  alone,  i.e.,  v  =  v(y),  the  wave  profile  is 
given  by 


(5.13)  f  =  q  3ech2  ^J^ 
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where 


K    = 


JJ^ 


0 


=  b(gd)' 


dy 


D 


^  n  [c  +  V  (y)]' 


and 


m  = 
o 


Tp   dzdy 
y 


D 


dzdy       p 


^  dy 
y  ^ 


-1  [c^+ v^(y)]' 


The  function  ^  is  defined  by 


o 


1 
~?  ' 

u 
o 


-1  <  y  <  0 


^y(O) 


u^{0) 
o^  ' 


^y(-l)    =    0 


ii{Q)    =   0  . 


Integration  gives 


and  so 


^y(y) 


r      d^ 


0 
%  =  bgd  J      [c^+  v^(y)]' 


y 


dt 


-1 


n   [C   +V   (?)] 

-1  '-  o   o 


2 

^  ^   dy 


Benjamin's  result  [13]  can  be  transformed  Into 


(5.14) 


2  ^1 


qK 


f  =  q  sech  ^  j^ 


where 
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0 
>2  r    dy 


K^  =  b(gd)-    ^ ^ 


and  -. 


Mq  :=  bgd  J   [v^(y)-  cf\[ 
-1  I  -1 


[v^(e)  -c] 


^-    ^  >  dy  . 


Here  c  denotes  the  actual  velocity  of  the  wave.   However,  as 
we  can  see  from  (4.51), 

c  +  c  =  0(1) 
o         ' 

which  means  that  (5. 13)  and  (5.14) 'give  the  same  first  order 
approximation.   Our  result  (5.13)  is  also  in -essential  agreement 
with  the  one  given  by  Ter-Krikorov  [12]. 

If  V  =  0,  (5.13)  reduces  to  the  familiar  result 

2  ""l  [3q 


(5.15)  f  =  q  sech^  ^   j^ 


It  is  easy  to  compare  this  with  the  results  corresponding  to 
^^[y)  =  C(y+1)  where  the  dimensionless  vortlcity  C  is  a  positive 
constant.   The  critical  speeds  for  constant  vorticity  are  given 

by 

0 

r      dy         1 


that  is 


^_i  [c^  +  C(y+i)]"^   sh 


Cq(Cq+  O  =  gh 


_  <±Jf+   4gh 
^o        5 
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or  In  terms  of  the  mean  speed  v  =  C/2: 


c   =  -V  ±  Igh  +  V 
o     o   ^^    o 


Also 


0  r  y  ,2 

m^  =  bgh  f     [c^+  C(y+1)]^)  f    ^ ^    dy 


-1  ^  -1  ^"o 


bgh 
and  we  have 


0 
K    ^  ,  2   r      dy 

""o       °  i^   [c^+C(y+l)] 


^ 


^Sh^o      /■       dy 


[c^+CCyi+i)]"^  i^   [c^+  C(y+i)]'^ 


IT  = ^ ?  '  -1  <  y  <  0 

""o   [c^ +C(yi+i)]^ 


by  using  a  mean  value  theorem.   If  c  >  0,    then 


"^o 


This  means  that  a  wave  of  given  amplitude  advancing  upstream 
on  a  stream  of  constant  vorticity  has  less  curvature,  is  broader, 
than  its  counterpart  on  a  stream  of  zero  vorticity.   If  c  ^  0, 
c  must  be  less  than  -C,  and  |  c^  |  >■  C;  so  for  this  case 


5^ 


3|c  1^ 
_K_  ^   ^'  o' 

"'o   [|c^|  -C(yi+i)]^ 


>   3 


which  shows  that  the  wave  moving  downstream  Is  sharper  than  Its 
counterpart. 

The  velocity  c  of  the  wave  Is  given  by 


_c   =   c^   +  1-^—  +   o   q) 

D        ^o 


Since 


K  -  b(gh) 


0 
2    r    ^dy 


.1  [c^  +  ay+i)]'^ 


_  b(3gh+  C^) 
3gh 


and 


D     %  -1  [Co+^fy+Dr 


we  have 


If   c^   >   0, 


b(2cQ  +  C) 
/gh 


c  -  c   +  q(^^  +  ^sh) 


4v2 

/ =^     _       q(gh+-^) 

-c   =  Jgh+  v^    -    v^    + 


o  o 


2high+  v^ 


o 
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If  c^  <  0, 

4V^ 


I —^       _  q(gh  +  -^) 

-c  =  -Jgh+  v^  -  v^ 


o     o       I — 

2h /gh  +  V 


o 


In  closing  we  remark  that  although  we  have  confined  our 
attention  to  solitary  waves,  we  can  easily  turn  the  analysis 
Into  one  on  cnoldal  waves  by  relaxing  the  conditions  (^-3^)  and 
(4.^8).   For  cnoldal  waves  we  find  the  equation 

NqTi!['(x)  -  N-^ti^(x)ti^^(x)  +N2ri-[(x)  , 

that  Is,  equation   (4.47)  but  not  with  the  same  coefficients. 
The  general  Integral  of  this  equation  satisfies 

This  leads  to  periodic  waves  of  the  form 

T]-j^  =  A  +B  cn(x,k) 

where  cn(x,k)  Is  Jacobi's  elliptic  function.   We  note  finally 
that  our  method  of  analysis  is  applicable  to  other  situations 
involving  long  rotational  and  irrotatlonal  waves  in  straight 
channels.   For  instance,  our  method  could  be  used  to  study 
solitary  waves  in  a  medium  of  two  layers,  with  different 
densities,  contained  in  a  channel  with  arbitrary  cross  section. 
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